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Abstract 

We consider elastic quark-quark scattering at high energy and fixed transferred momen- 
tum. Performing factorization of soft gluon exchanges into Wilson lines vacuum expectation 
values and studying their properties, we find that the asymptotics of the scattering ampli- 
tude is controlled by the renormalization properties of the so called cross singularities of 
Wilson loops. Using this fact, we evaluate the scattering amplitude and show that its asymp- 
totics is determined by the properties of the 2x2 matrix of anomalous dimensions which 
appears after one renormalizes the cross singularities of Wilson loops. A generalization to 
the case of quark-antiquark and gluon-gluon elastic scattering is discussed. 
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1. Introduction 



Investigation of the hadron-hadron interaction at high energy s and fixed transferred momentum t is one of the 
longstanding problems in QCD ^. The first calculations of the elastic near forward parton-parton scattering 
amplitudes to the lowest orders of perturbation theory (PT) in the leading logarithmic approximation revealed 
novel features of nonabelian gauge theories It was found in QCD that gluons, in contrast with photons 
in QED Q, are reggeized and their contribution to the scattering amplitude has the asymptotics ~ s"^*-' with 
the Regge trajectory a{t) known to the lowest order of PT. Contributions of individual Feynman diagrams to 
the scattering amplitude obey a remarkable property. In the center of mass frame of the incoming particles 
the s— dependence of diagrams comes from integration over longitudinal components of the exchanged gluon 
momenta whereas the dependence is given by two-dimensional integrals over transverse components of gluon 
momenta. It was proposed ^ to explore this property in order to describe the asymptotics of the scattering 
amplitudes within the framework of some effective two-dimensional field theory. Lipatov succeeded to resum 
the so called itt— corrections using explicit expressions for the S*— matrix of exactly solvable two-dimensional 
field theories ^. Most recently, another two-dimensional field theory was proposed to describe high-energy 
behavior in QCD [Q. Although this theory reproduces the known expressions for the amplitude to the lowest 
orders of PT in the leading logarithmic approximation, it is very hard to solve the model exactly to find the 
asymptotics of the scattering amplitude to all orders of PT. 

The asymptotics of the parton-parton scattering amplitude was found in the leading logarithmic approxi- 
mation using the "evolution equation" approach |0, |[ ||]. In this approach, one uses the results of one-loop 
calculations of the scattering amplitude in order to resum the leading logarithmic corrections to all orders of 
PT. However, in spite of considerable progress, we don't have a regular way to resum and control nonleading 
logarithmic corrections to the scattering amplitude. In the present paper we propose a new approach which 
enables us to take into account systematically these nonleading corrections. 



2. Factorization of soft gluons 

Let us consider near forward hadron-hadron scattering with large energy s and fixed transferred momenta t 
which is chosen to be much larger than the QCD scale Aqcd- In this kinematics, interaction between hadrons 
occurs at short distances and we may apply the parton model to relate the asymptotics of the process to 
properties of the parton-parton scattering. We start by considering in detail a near forward elastic quark-quark 
scattering and generalize consideration to quark-antiquark and gluon-gluon scattering at the end of the paper. 

Let pi,p2 and p'i,p'2 be the momenta of incoming and outgoing scattered quarks, respectively. We choose 
all these momenta to be on-shell Pi=P2= p'l = p'2 = ™^ with m being quark mass and define the kinematic 
invariants 

s = (Pl+P2)^ t={pi~p[f, u={pi-p'2f 

which are related as s + t-\-u ~ Am? . The amplitude of elastic quark-quark scattering contains color indices 
of both incoming and outgoing quarks. It depends on these invariants as well as on the quark mass 

m and infrared cutoff A (a fictitious gluon mass, for instance) one has to introduce to regularize IR divergences. 
We will study the asymptotic behavior of the amplitude of the near forward elastic quark-quark scattering in 
the following region of parameters: 

s, TO^ > > > Aqcd- 

The amplitude T*^"' is a function of the ratios of the kinematic invariants. One could propose a lot of such 
ratios but, as we will show below, only one possibility is realized:^ 

Moreover, the s— dependence of the amplitude comes through the dependence on the angle 7 between quark 
4— velocities Vi = pi/m and V2 = Pij™, defined in Minkowski space-time as 

(wi • V2) — cosh7, s — 2rr?(\ -I- cosh 7) ~ Am"^ cosh^(7/2) . 

"'^This dependence is confirmed by calculations |^ of the scattering amplitude to the lowest orders of PT 
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Let us go to the center of mass frame of the incoming quarks. In this frame, quark momenta have the following 
light-cone components: 

Pi = ipt,Pi,Pi) = ^(e^/^e-^/^0), P2 = ipt,P2,P2) = ^(e'^/^e^/^0). 

In the limit s 3> the angle between quark velocities becomes large 7 = log ^ 3> 1 and both quarks move 
close to the "+" and "— " light-cone directions. The total momentum transfer is q = pi —p'l, and in the center of 
mass frame its components are (j+ = —q^ = -^=-^===5 — 0{t/y/s) and q^ = tu/{s — 4m'^) = 0{t). Thus, in the 

limit s 3> ~t, we can neglect the longitudinal components of the transferred momentum and put q ~ (0^, 0^, q) 
with q^ = t — ~q^. 

Note, that the transferred momentum is much smaller than the energies of incoming quarks. This suggests 
that the incoming quarks interact with each other by exchanging soft gluons in the t— channel with the total 
momentum q. This interaction gives rise to infrared divergences of the scattering amplitude. It is a general 
feature of soft gluons Q that their contribution to the scattering amplitude is factorized into a universal factor 
given by a Wilson line vacuum expectation value. The origin of this property is the following. Interacting with 
soft gluons quark behaves as a relativistic charged classical particle, and the only effect of its interaction with 
soft gluons is the appearance of an additional phase in the quark wave function. This phase, the so called 
eikonal phase, is defined as a Wilson line V exp(i dx ■ A{x)) evaluated along the classical trajectory C of the 
quark. Applying the same approximation to the near forward elastic quark-quark scattering we get two Wilson 
lines corresponding to two incoming quarks. Note, that although quarks change their velocities after scattering, 
we can neglect this difference in the limit —t ^ m^. It means that calculating the eikonal phases we may 
consider incoming quarks as classical particles moving from —00 to +00 with velocities vi and V2- We combine 
two eikonal phases of both incoming quarks and get the expression for the scattering amplitude in the form 

Tff' (^-^, = sinh7 J d^ze-'^-'^{0\Twi\0)W'yiz)\0) , t = -f (2) 

Here, Wilson lines W+ and W- are evaluated along infinite lines in the direction of the quark velocities vi and 
V2, respectively: 



W+(0) ^Vexp (^i J davi- A{via)^ , W-{z) = "Pexp (^i J 



dp V2 ■ A{V2P + z) , (3) 



and the integration paths are separated by the impact vector z — (0^ ,0 , z) in the transverse direction as 
shown in fig.l. 




Fig. 1: Wilson lines are evaluated along classical trajectories of quarks with velocities vi 
and V2 separated by the impact vector z = (0"*", 0~, z) in the transverse direction. 

Integration over two-dimensional impact vector z was introduced to ensure the total transferred momentum 
to be equal to q = (0+, 0^, g). Indeed, by expanding the Wilson lines in powers of gauge fields and using the 
momentum representation for these fields one gets that expression (^, firstly, reproduces the eikonal approxi- 
mation for interaction vertices of incoming quarks with soft gluons and, secondly, the total momentum of gluons 
ktot. propagating in the channel is restricted by the following 5— functions: 

S{ktot. ■ vi)5{ktot. ■ V2)5{ktot. - q) = -A. — <5(fcj+j )6{kiot )S{ktot. - q) = -t\ — 5{ktot. - q) 

smh7 ■ ■ smh7 
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where the first (5— function comes from W-f, the second one from W- and the last one from integration over 
z. To compensate the additional factor in the r.h.s. of this relation, the same factor was introduced in the 
expression (||). 



3. Properties of Wilson loops in QCD 

The calculation of vacuum expectation values of Wilson lines entering into the expression (^) for the scattering 
amplitude is a hard problem in PT.0Let us consider first a special case of the elastic electron-electron scattering 
in QED. One may use the expression (||) for the scattering amplitude in this case, but the calculation of the 
Wilson lines is simplified because for abelian gauge group a path-ordered exponent coincides with an ordinary 
exponent. Moreover, since photons don't interact with each other, the calculation of Wilson lines in QED is 
reduced to the integration of a free photon propagator along the path of fig.l and the result of calculation is 

(0|TM^+(0)M^_(z)|0) =exp(^-y dx'' J dy"" D^,{x - y)^ = (a^z 2/4)*^ (4) 

where A is a fictitious photon mass. After substitution of this relation into and integration over impact 
vector we get the expression for the scattering amplitude in QED [Q 

rf^ ±\ = ■ /_A-i-"'^ r(i-i-^|icoth7) 
W'xy t UV r(l-z|icoth7) 

Of course, there are no reasons to expect that calculation of the Wilson lines in QCD will be simple, but there 
are some features of the scattering amplitude in QED which have their analogs in QCD. In particular, the 
scattering amplitude obeys the condition ([^) and its s— dependence comes from the dependence of Wilson lines 
on the angle 7 between electron velocities. In the limit s/rri^ — s- 00 or, equivalently 7 — > cx), the expectation 
value (^) doesn't depend on s and, as a consequence, the scattering amplitud has the asymptotics T/TBom 
which implies that photon is not reggeized in QED. Here, Tsom = the expression for the amplitude in 

the Born approximation. 

The asymptotics of the scattering amplitude (||) in QCD is defined by properties of the line function equal 
to the vacuum expectation value of the product of two Wilson lines: 

Wff' = {0\TWl\0)Wi' {z)\0) (5) 

It depends on the color indices of incoming and outgoing quarks and is gauge invariant. According to the 
definitions (^) and (^, the line function W depends, in general, on the quark velocities vi and V2, the impact 
vector z and infrared cutoff A which is introduced into Feynman integrals to regularize IR divergences. Trying 
to form dimensionless scalar invariants formed by these variables we find that, due to identities vf ^ V2 = 1 
and (viz) = {v2z) = 0, there is only one possibility 

= A^^') ^ W;^'(7, A^^^) (6) 

After substitution into (^), this relation implies the functional dependence (|]) of the scattering amplitude. 
Thus, the s— dependence of the scattering amplitude comes from the dependence of the line function on the 
angle between quark trajectories in fig.l while its i— dependence is related to the dependence of the line function 
on the impact vector z. 

To understand the z— dependence of line function, let us consider as an example the one-loop calculation of 
W in the Feynman gauge. Using the definitions (^) and (||), we get 

2 /"OC /"OO / \ 

W^.ioop = I^I+it''^n-^r{D/2~l)X'-^ da dl3— r,'"! -.2 , ...n/.-i (7) 



47r^/2 i_oo i-00 [-(i'ia-V2/3)2+z2+tO]^/2 

where a direct product of the gauge generators defined in the fundamental representation takes care of the color 
indices of the quarks. In this expression, gluon is attached to both Wilson lines at points via and V2(3 + z 



^The effective two-dimensional field theory proposed in enables us to calculate this object only to the 
lowest orders of PT 
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and we integrate the gluon propagator ^2 -D^i/(f la — W2/3 — z) over positions of these points. To regularize 
IR divergences we introduced the dimensional regularization with D = 4 + 2e, {e > 0) and A being the IR 
renormalization parameter. There is another contribution to Wi-ioop corresponding to the case when gluon is 
attached by both ends to one of the Wilson lines. A careful treatment shows that this contribution vanishes. 
Performing integration over parameters a and (3 we get 

Wi-ioop{l, A^z") e) — {-n:cot\Y-i)T{e){^\^z\' (8) 

TT 

The integral over a and (3 in (|^) has an infrared divergence coming from large a and (3. In the dimensional 
regularization, this divergence appears in Wi-ioop as a pole in (D — 4) with the renormalization parameter A 
having a sense of an IR cutoff. This result may seem strange because, after substitution of the one-loop expression 
for W into (0), we should reproduce the expression for the scattering amplitude in the Born approximation 

■ 2 

TBorn = ^;^(^° ® corresponding to one-gluon exchange in the i— channel. Indeed, performing a Fourier 
transformation of Wi-ioop in (j^ we find that IR divergence disappears, and the one-loop expression TBom for 
the scattering amplitude is reproduced. 

The approach we propose for the calculation of Wilson lines vacuum expectation value is based on the 
following observation. The one-loop expression (||) for the line function is divergent for z = 0. This divergence 
has an ultraviolet origin because it comes from integration over small a and /3 in that is from gluons 
propagating at short distances (wia — U2/3) between quark trajectories. One should notice that the ultraviolet 
(UV) divergences of Wilson lines at z = have nothing to do with the "conventional" ultraviolet singularities 
in QCD. For z = 0, the integration paths of Wilson lines of fig.l cross each other at point and UV divergences 
appear when one integrates gluon propagators along integration path at the vicinity of the cross point. It is not 
for the first time when one encounters these very specific divergences of Wilson lines. The same divergences, 
the so called cross divergences, were observed more than 10 years ago [ p^ when attempts have been made to 
reformulate the dynamics of gauge fields in terms of string operators. 

We note that the cross divergence has not appeared in Wi-ioop for nonzero z because, as it enters into the 
integral (0), nonzero z ^ regularizes gluon propagator at short distances. Thus, z ^ has a meaning of an UV cutoff 
for the one-loop line function W and the z— dependence of the original line function VF(7, A^z ^) (or, equivalently, 
the i— dependence of the scattering amplitude) is in one-to-one correspondence with the dependence of the same 
Wilson line but with z = on the UV cutoff which we have to introduce to regularize the cross singularities. 
The general reason for this property to be valid to all orders of PT is the following. As follows from (^, the line 
function depends only on two dimensionfuU arguments: z^ and A^ while quark velocities are dimensionless. 
It implies that Feynman integrals over gluon momenta which one gets by calculating vacuum averaged Wilson 
lines in the perturbation theory contain only two momentum scales: A^ and 1/z^. The fact that the first scale 
cuts gluon momenta from below implies that the second scale cuts gluon momenta from above. For z = 0, this 
second scale goes to infinity and line function becomes UV divergent. For z^ 0, the line function is finite 
but the cross singularities manifest themselves in the singular dependence of W on z as z ^ ^ 0. The classical 
trajectories of quarks in fig. 1 are shifted by vector z in the transverse direction, and we may consider this shift 
as one of possible ways to regularize the cross singularities which line function W will have for z = 0. This 
simple observation suggests the following way for calculation of the z— dependence of the Wilson lines: first, put 
z = in the definition of line function and introduce the regularization of cross singularities of W] second, 
renormalize cross singularities and identify UV cutoff with impact vector as /i^ = 1/z^. 

Thus, the z— dependence of the line function W and, as consequence, the asymptotics of the scattering 
amplitude are controlled by the renormalization properties of the cross singularities of Wilson loops. 



4. Renormalization of Wilson loops 



The line function W^j has color indices of incoming and outgoing quarks. Since quarks are defined in the 



fundamental representation of the gauge group SU{N) one can decompose W into two invariant tensors darSjji 
and Sij'Sji' : 



Projections of W onto these tensors define two gauge invariant Wilson loops Wi"'^ = ^ij'^ji'^i/ W^'^^ = 
Sii'Sjj'W^j'' which are evaluated along the integration paths of fig. 2 (a) and (b), respectively, closed at 
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infinity. According to the approach formulated at the end of the previous section, we put z = and consider 
renormaUzation properties of these Wilson loops having a cross point. 




(a) (b) 

Fig. 2: Integration paths entering into the definition of the "physical" Wilson loops VF^*'"'' 
and W^''\ respectively. Dotted line indicates the way in which these paths are closed at 
infinity. 

It turns out that the renormalization properties of Wilson loops with a cross point depend on the repre- 
sentation of the gauge group in which Wilson loop is defined. In the case of the quark-quark scattering, the 
Wilson loops are defined in the quark representation but considering a near forward gluon-gluon scattering in an 
analogous manner we will get the same expressions for the scattering amplitude but with Wilson lines defined 
at the adjoint representation of the gauge group. It is exactly the way in which we will find the difference 
between the asymptotics of the amplitude of quark-quark and gluon-gluon scattering. 

For Wilson loops Wi""^ and W^''^ defined in the fundamental representation of the SU{N) group the renor- 
malization properties can be formulated as follows 0. The Wilson loop W^"^ of fig.2(a) is mixed under 
renormalization with the Wilson loop W^j"-* of fig. 3(a), and the Wilson loop W^''^ of fig. 2(b) is mixed with the 
Wilson loop I^f^ of fig.3(b). 




(a) (b) 

Fig. 3: Wilson loops VFj"'' and Vl^2^'' which are mixed under renormalization with the "phys- 
ical" Wilson loops of Fig. 2(a) and (b) respectively. 

Both integration paths in fig. 3(a) and 3(b) are different from the original paths of fig. 2(a) and 2(b) only at the 
vicinity of the cross point. It is natural to introduce the following doublets: 

--CI) 

which don't mix with each other under renormalization.^ Then both doublets obey the same renormalization 
group equation (lO| 

+ /3(5)^) W^(7, AVm') - -Tcrossh, g)W{-,, AVm') (10) 



^Notice that the Wilson loops W^"^ and W^'''' were introduced into consideration by the renormalization 
procedure and in contract with W^'^^ and Wi''^ they don't have any physical meaning 
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where ^ is the renormaUzation parameter and A is an IR cutoff. A new object appears in the r.h.s. of this 
equation, a 2 x 2 matrix of anomalous dimensions Tcrossi'J, g)- This matrix is gauge invariant and it depends 
only on the coupling constant and the angle between quark velocities at the cross point. Following [lO|] we 
calculate the one- loop expression for this matrix in the fundamental representation of the SU{N) gauge group 

r,rossh,9)^-m + 0{al), (11) 

TT 

where 2x2 matrix r(7) depends only on the angle 7 and is given by 

Y( ) — f ^ifcoth7 i7rcoth7 

\^ — 7C0th7 + 1 + i7rcoth7 Af(7 coth7 — 1) — ^ coth7 

The RG equation ( [lo| ) implies that the s— dependence of Wilson loops is defined by the properties of this matrix. 
The one-loop expression (^l]) is valid for an arbitrary values of 7. In the limit s ^ m^, the matrix r(7) has the 
following asymptotics: 

ITT 

-l0g;^ + l+Z^ 7V(l0g;^-l) - 



r(7) = ^ ^ 

V - log + -1- -r ly ^1^5 - ---J - N 

where we neglected terms vanishing as s/m^ -^00. To check the RG equation ( p^ ) we perform the one-loop 
calculation of Wilson loops 

VKi'"^(7,AV/i') = A^fl- — CF•2^7^coth(7)log^) 

V TT A / 

Wt\j,xy^,^) = iV2(l-^C^.2(7Coth(7)-l)log^) 
iyf)(7,AV/i^) - iV2(l-^-01og^) 

W^f^(7,AVM') = ^fl- — CF•2^7^coth(7)log^) (12) 

with Cf — ^2N^ ^'^^ verify that they do satisfy the RG equation (p^). Note, that these Wilson loops have 
different normalization at the lowest order. Moreover, at /i = A the one-loop corrections disappear 

M^('')(7,AVm^)U=a = (^,), 1^(^)(7,AVa^^)U=a = C;^') (13) 

which means that for UV and IR cutoff for giuon momenta equal to each other there is no phase space for 
gluons. 



5. Solving the RG equation 

By solving the RG equation ( pO[ ) with the boundary conditions (|l^) we could find the "physical" Wilson loops 
Wi"^^ and Wi^\ Different Wilson loops of figs. 2 and 3 satisfy the same equation (|l^) and the reason why we 
get four different expressions for them lies in the different boundary conditions (|l^) we impose on the solutions 
of (p^). The general solution of the RG equation ([l^) with the boundary conditions ( p^ ) for the doublet ly*^") 
is given by T— ordered exponent 

IF('^)(7, AVm') = Texp £ ^r_,3(7, 5(r))) ) , (14) 

and for the second doublet we get analogous expression 

PF('')(7, XV n') = Texp J' ^T^rossh, 5(r))) (^^') . (15) 

The matrices Tcrossi'y, <?(''■)) don't commute with each other for different values of r and it makes very difficult to 
evaluate the T— exponents]^ To find the scattering amplitude we extract the upper components of the doublets 

^We may follow another route and eliminate one of the components of the doublet W from the system of 
equation (|lO| ) to get the equation for the second component. The resulting second order differential equation is 
known as the Riccatti equation 
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( |l4[ ) and (|l5|) and get two "physical" Wilson loops Wi""^ and Wi''\ evaluate the Wilson lines VF/ using (^), 



identify the renormalization parameter as = 1/z^ and substitute W^^-' into (|^) 



where An and A12 are elements of the 2x2 matrix ^4.(7, z^A^) = Texp f — J^^^^ ^rj,ross(7j j ■ 

The simplification occurs when we substitute the one-loop expression (|ll|) for the matrix Tcross into (|lj) and 
(|l5|). One-loop matrices rcross(7i commute with each other for different r and we can omit T— ordering 
in the definition of the matrix A. After diagonalization of the one-loop matrix r(7) by a proper unitary 
transformation we get 

^(7,Az) = ^--^ exp(^-r_y^ -^j+Y—Y7'^^[-^ 
Here, T± are the eigenvalues of the matrix r(7) which satisfy the following characteristic equation 

r± - r± ^iV(7 coth 7 - 1) - ^ coth 7^ + tt^ coth^ 7 

Finally, we find the following expression for the scattering amplitude ( p^ 

Ti/' = S,,,S,,, t("' + 4^,^ t(«) (17) 

where two invariant amplitudes corresponding to the exchange in the i— channel by states with the quantum 
numbers of vacuum and gluon are given by 



„/n-i sinh7 r_i_r^ „ , „/o-, cosh7 r-i_T_i_ — r_r_ 

t(o) = — _j:_^_(t+ - T_), = 2zn ^ ^ r+ - r_ ^^^^ 

Here the notation was introduced for the Fourier transforms 

T,^ A /,.,e--^^expf-r, r^'!t!^]^± (d^ze---^C':^'lf,fr 



-2r±/A, 



To the lowest order of PT we have T± = 2as and after substitution of these values into we recover one-loop 
expression for the scattering amplitude: T(°) = and T^*^) = To get the expression for T± in the leading 

log approximation we freeze the argument of the coupling constant and after integration get 

T± ~ 2a, exp (-^r± log ^] . (19) 



2tt ^ ° X 



Let us consider the properties of the obtained expressions (|T^ in the limit of high energies s 3> ■ Solving the 
characteristic equation we find that in the large s limit both eigenvalues have positive real part (Re T± > 0). 
Moreover, one of the eigenvalues of matrix F is much larger than the second one: 



2^ - 1 , -1 s . ^9 s 



N3 



log-i4+o(log^^4) 



Together with ( p^ ) this property of the matrix Tcross implies that the amplitudes of singlet and octet exchanges 
have the following high-energy behavior 

T^°^ cx r_r+ - F_r_ , t'^^ cx t+ - — t_ . 

r+ 

The functions r+ and T- have a Regge like behavior and for F+ ^ F_ we have r+ <C T_. Moreover, the 
contribution of T_ to the amplitude of the octet exchange is suppressed by the factor l/F^ = ©(log"^ s/m^) 
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as compared to that to the amphtude of the singlet exchange. Thus, the high-energy asymptotic behavior of 
the scattering amphtude ( p^ ) is dominated by the contribution of r_ to the singlet exchange amplitude T^'^'' . 

To get the expression for the scattering amplitude ( [l7| ) in the leading logarithms in s and t we substitute 
r+ = TV log s/m^ and r_ = into (|l8|) and (p^). The result has the standard reggeized form T = 
TBo™(s/m2)«(*) with the Regge trajectory a{t) = -^iVlog^ = -^N J -l^r^. 

6. Conclusions 

In this paper we considered the elastic quark-quark scattering at high energy s and fixed transferred momentum 
t. Performing factorization of soft giuon emissions, we found the expression for the scattering amplitude (H) as 
a Fourier transformed vacuum expectation value of a product of two Wilson lines defined in the quark represen- 
tation of the SU (N) gauge group and evaluated along the classical trajectories of quarks with velocities vi and 
V2- We calculated these Wilson lines using the fact that in the limit of zero impact vector z, when quark trajec- 
tories intersect each other, Wilson lines get very specific UV cross divergences. There is a close relation between 
the renormalization properties of cross divergences and the z— dependence of Wilson lines or, equivalently, the 

dependence of the scattering amplitude. Using this relation we found the scattering amplitude ( p^ which 
depends on the 2x2 matrix of anomalous dimensions Tcrossili g) of Wilson loops. Performing the one-loop 
calculation of this matrix we found that these are the properties of the matrix Tcrossij, s) which determine the 
asymptotics of the scattering amplitude (0). Being a function of the angle 7 between quark velocities, the 
matrix Tcross{l, g) depends on the energy s and one has a regular way |]l0| for its calculation in higher orders of 
PT. The important question arises: whether higher order corrections will change the asymptotics of the one-loop 
matrix ( |ll|) and the properties of the eigenvalues T± which were important for us to get the asymptotics in 
(|l7|). This will be the subject of a separate paper. 

We may generalize the consideration to study the high-energy elastic quark-antiquark and giuon-gluon 
scattering. In both cases, the scattering amplitude is given in a close analogy with (||) by the vacuum expectation 
value of two Wilson lines coming as eikonal phases of incoming particles. It is important to recognize that Wilson 
lines carry the same representation at which these particles are defined. For quark-antiquark scattering, the 
Wilson line coming from antiquark is evaluated along its velocity V2 , say, with gauge generators in the antiquark 
representation. We notice that in this situation an antiquark can be treated as a quark moving backwards in 
time with velocity —V2 and all previous results for the quark-quark scattering amplitude can be applied provided 
that we replace one of the quark velocities as V2 — * ~V2- In terms of the angle 7 between quark velocities this 
transformation corresponds to the replacement 7 — > jtt — 7 in Tcrossil, g)- For gluon-gluon scattering both 
Wilson lines are defined in the adjoint representation of the SU{N) gauge group. Following the analysis of 
Sect. 3, we relate the asymptotics of the gluon-gluon scattering amplitude to the renormalization properties 
of Wilson loops of fig. 2. As it was stressed before, these properties are different for the fundamental (quark) 
and the adjoint (gluon) representations of the SU{N) gauge group. One may find a novel feature of cross 
singularities in the adjoint representation by considering the one-loop expression (H) for the Wilson lines in 
the simplest case of the SU{2) gauge group. Simphfying the direct product of the gauge generators of the 
SU{2) group we get the following color fiow: tf,^tj,j = \5iji5iij — jSi'iSjij in the fundamental representation and 
^i'i^fj — ^ij'^i'j ~ SijSi'j' in the adjoint representation. Different color flows in these two representations imply 
that Wilson loop of fig. 2(a) defined in the fundamental representation of SU(2) is mixed under renormalization 
with the Wilson loop of fig. 3 (a) while in the adjoint representation there is an additional Wilson loop of fig. 4 
which is mixed with both of them. 
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Fig. 4: Wilson loop which is mixed under renormalization with the Wilson loops of Fig. 2(a) 
and 3(a) defined at the adjoint representation of the SU{2) group. 

Coming back to the SU {N) group, we note that the renormalization of cross singularities was performed |l0| 
only in the fundamental representation and the generalization to other representations is an open question. 

There have been a number of attempts in the 80's to describe the confining dynamics of QCD in terms 
of Wilson loops. One of the problems which has arisen on this way was jll], |l^ that Wilson loop possesses 
extra UV singularities when the integration path has end point, contains cusp or crosses itself at some point. 
The renormalization of these singularities introduces into consideration new anomalous dimensions of Wilson 
loops: Tendig), ^cuspijjg) and Tcrossil, g), respectively. It turns out that these anomalous dimensions 
which appear as undesired features of Wilson loops are of the most importance in perturbative QCD when one 
studies effects of soft ghions. Namely, the infrared behavior of quark and giuon propagators is controlled by 
Tend- The asymptotics of the Sudakov form factor, the structure functions of deep inelastic scattering for x —^ 1, 
large perturbative corrections to the Drell-Yan cross section are determined by Tcusp, the same function was 
introduced in the heavy quark effective field theory as the velocity dependent anomalous dimension. Finally, as 
it was shown in this paper, it is T cross ^ which governs the high energy behavior of the parton-parton scattering 
amplitudes. 
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